Equation of motion of an interstellar Bussard ramjet with radiation and mass losses 
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An interstellar Bussard ramjet is a spaceship using the protons of the interstellar medium in a 
fusion engine to produce thrust. In recent papers, it was shown that the relativistic equation of 
motion of an ideal ramjet and of a ramjet with radiation loss are analytical. When a mass loss 
appears, the limit speed of the ramjet is more strongly reduced. But, the parametric equations, in 
. . . terms of the ramjet's speed, for the position of the ramjet in the inertial frame of the interstellar 

, medium, the time in this frame, and the proper time indicated by the clocks on board the spaceship, 

■ can still be obtained in an analytical form. The non-relativistic motion and the motion near the 

' limit speed are studied. 
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I. INTRODUCTION 



. . During its motion through the space, an interstellar fusion ramjet collects interstellar ions (mostly protons) with a 

' magnetic scoop (or ramscoop) to supply a fusion reactor able to fuse protons to obtain helium [U, 0, y, [3| ■ Accelerated 

II reaction products are exhausted out of the spacecraft's rear to produce thrust. An ideal ramjet could reach a velocity 

^ very close to the speed of light c. If some energy extracted from the interstellar medium is lost in form of thermal 



radiation, the ramjet speed is limited to a value below c f5|. In this two cases, we have shown that analytical formulas 
can be obtained for the position of the ramjet in the inertial frame of the interstellar medium, the time in this frame, 
and the proper time indicated by the clocks on board the spaceship 0| • These parametric equations are given in 
terms of the ramjet's speed. 

Moreover, it is natural to assume that a fraction of the collected interstellar gas can be lost during the work of the 
engine. The limit speed is then more strongly reduced [sl and the equation of motion further more complicated. In 
this paper, we show that analytical formulas can also be obtained for the position, the time and the proper time of 
^ i' a ramjet with radiation and mass losses. The non-relativistic motion and the motion near the limit speed are also 
studied. 

^ ' The situation presented in this paper is more realistic than the case of a Bussard ramjet with a perfect fusion engine 
, 0- Nevertheless, it is not very probable that the kind of ramjet considered here could be ever built, even in the far 
' future. Consequently, this work can be considered as an advanced exercise (from the point of view of calculus) in 
^\ , special relativity. The framework of interstellar space travel could be very attractive for undergraduate students. The 
■ basic equations are simple outcomes of momentum and energy conservation. Even if the solutions demand very heavy 
C5 calculations, they have interesting properties to explore. However, let us note that the ramjet concept is potentially 
too valuable to be simply discarded despite its tremendous technical difficulties. Some researchers have suggested 
T-H alternatives to the initially proposed proton-proton fusion ramjet [1,01: the catalytic ramjet (the use of catalyzed 
C . fusion reaction with a high rate), the RAIR (the use of nuclear fuel carried by the ship), etc. 



II. GENERAL EQUATIONS 



' In the following, all calculations are performed in the frame of the interstellar medium, considered as an inertial 
frame. A Bussard ramjet of constant mass M moves at speed v = (3c through this medium, which contains protons at 
rest with a mass density p. The effective intake area of the ramscoop is denoted A. A fraction e of the absorbed mass 
is converted into useful kinetic energy in the hydrogen fusion reactor, a fraction A is dissipated in form of thermal 
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radiation and a fraction k is lost in the interstellar medium (parameters e, A and k have the same meaning than in 
Ref. [^). This means that, if a mass dm of protons at rest is scooped up from the interstellar medium, only an energy 
e dm (? is converted into ordered motion of the exhausted material. 

In the frame of the spaceship, we can assume that the thermal loss is isotropic and that the radiation carries no 
momentum. For the ramjet, the energy absorbed due to the mass dm is dE = dm^ ^ where 7 = 1/ ^/l — fp , and the 
engine dissipates an energy equal to XdE. Lorentz transformations [1, Q imply that an energy ^XdE — Xdmj^c^ 
and a momentum 7 /3 XdE/c = Xdm^'^/Sc are lost in the interstellar medium frame. 

If a fraction k of the collected interstellar gas is lost during the work of the engine, a useless mass k dm is dropped 
out with no velocity in the frame of the ramjet. So, this mass has an energy ndm^c^ and a momentum ndm^fic in 
the interstellar medium. 

During the time interval dt^ a mass dm of protons, at rest, is scooped-up. The ramjet speed is then increased by 
the quantity cd(3^ thanks to the ejection of a mass (1 — a) dm of helium with a speed wc, where 

a = e + A + K. (1) 

Another interesting quantity is the fraction of matter, 

uj = e + X, (2) 

which is converted into pure energy by the engine. 
The conservation of momentum implies that 0, Q 

M 7(/3) 13c = M "f{f3 + dp) {(3 + d/3)c + A dm 7(/3)^/3 c + k dm 7(/3)/3 c + (1 - a) dm 7(10) wc, (3) 

where 7(x) = l/yl^^a? and dj3 > 0. The conservation of energy leads to 

M 7(/3) + dm = Af 7(/3 + d/3) + Xdm-/{Pfc^ + Kdm-/{f3)c^ + {1 - a) dm -f{w) . (4) 

The collected mass is a function of the ramjet speed and is given by 

dm = A p f3c dt. (5) 



III. ACCELERATION 



Taking into account equations it is possible to compute the acceleration Lp = dv/dt oi the ramjet, measured 

in the rest frame of the interstellar medium. With this aim, it is natural to define characteristic acceleration 1^9*, time 
and length x<t by the following relations 

M ip^ ip^ 

It is also useful to introduce the following notations 

u' — u{2 — u), u = 1 — u, (7) 

with 1 — u' = V?. 

It is worth noting that u is only 0.0071 for the most energetic known fusion reaction [1]. So, for all fusion reactions, 
we have « 2ijj and a) « 1. If the particles collected by the ramjet are an ideal mixing of matter and antimatter, 
the mass reaction can be totally converted in pure energy {uj = 1 and k = 0) [6^. 

Using the relation 

7(/3 + d/3) =7(/3)+7(/3)3/3d/3, (8) 

the elimination of the reduced speed w of the exhausted reaction mass by the relation 7(w)^(l — w"^) — 1 gives a 
second degree equation in ip, whose physical solution is (in the following equations, the simplified notation 7 means 
more precisely 7(/3)) 



^^-211 ( V72 - 1 + F(7) - Vl^) 



with (9) 



i^(7) = -A'72 - 2kA7 + k2 (10) 



FIG. 1: Limit reduced speed f3i of tiie ramjet as a function of the parameters e, A and k. Surfaces colored in dark grey to light 
grey correspond respectively to /3; = 0.2, 0.4, 0.6, 0.8. Left: Arbitrary values of e, A and k. Right: Values of w < 0.0071 with 
the vertical white plane marking the constraint uj = 0.0071. 



This equation gives the acceleration of the ramjet as a function of its speed in the inertial frame of the interstellar 
medium. If a = (e = A = k = 0), the acceleration vanishes as there is no input of energy into the engine. Since 
</? = if /3 = (7 = 1), an initial boosting is necessary for the ramjet. This is due to the fact that the reaction 
mass reaches the reactor thanks to the speed of the ramjet. In theory, a very small speed is sufficient to start the 
ramjet. In practice, a fusion reactor could probably not operate correctly without a sufficient intake. Nevertheless, 
the ramjet could accelerate with an initial speed as low as 10 km/s Such a speed could be reached with usual 
chemical rockets or by future nuclear rockets. As expected, the acceleration also vanishes for f3 — I ("f — 00) since no 
object can move at the speed of light. 

But this limit speed is never reached for a non-ideal ramjet. It is clear from Eq. Q that ip — when -^(7) = 0. 
F^j) is a quadratic function in 7 with two roots — F; and 7; such that 

F(7) = A'(rz+7)(7i-7) with (11) 



Vk2 + X'a' - Ak , + X'a' + \k , , 

11 = y and = . (12) 



With a little algebra, it can be shown that 



1< 7i < Ti, (13) 



taking into account the condition a < 1 to obtain a non- vanishing exhausted mass. For relevant values of the speed 
(7 ^ l)j one can see that: F{'ji) = 0, ^(7) > for 7 g [1, 7;] and ^(7) < for 7 > 7/. So, 7/ is a limit value: if 7 = 7;, 
the ramjet acceleration vanishes (</? = 0); if 7 < 7;, the ramjet speed increases ((^ > 0); if 7 > 7;, the ramjet speed 

decreases {(p < 0). Thus, the ramjet speed tends toward the limit reduced speed /?; — ^1 — In the following, we 

will only consider the realistic case (3 < f3i or "f < ji. Contrary to the case of an ideal ramjet, the velocity of a ramjet 
with radiation and/or mass losses cannot be arbitrarily close to the speed of light [1] (see Figs. [T] and [2]). 



Even with small values of the parameters A and k, the maximum speed can differ significantly from 1. This means 
that the slowing down of time on board the spaceship can become not large enough to allow interstellar travels in a 
period of time bearable for human beings. For values of A + k 3> e, the motion of the ramjet becomes non-relativistic. 

In our previous papers [1, all equations were given as function of the reduced speed (3. We find here more 
convenient to present the parametric equations of motion of the ramjet as a function of its reduced speed through the 
kinematical factor 7 = 7(/3). We will assume that, at a time t = in the inertial frame, the position of the ramjet is 
X — and its reduced speed is (3o ^ 0. Moreover, the clocks on board the ramjet indicate a proper time r = 0. 



IV. TIME 



Since 



dt dt c 

the solution of Eq. Q is given by the foUowing integral 



^ T P-TT = 7 V 7 - (14) 



^ = r , ^. (15) 



t* Jo J JO (72 - 1) (v/7^ - 1 + F(7) - V7^~l) 
We give here the main steps of the procedure to solve this integral: 



• To multiply the numerator and the denominator of the fraction by the quantity (y^72 — 1 + ^"(7) + 1/72 — 1) 
to obtain the function (72 — 1)^(7) at the denominator; 

• Using the relation (jlip . to write the fraction obtained as a sum of simpler fractions; 

• To integrate each of these new fractions and to simplify the result. 

The calculation is very heavy, but an analytical form can be found. To get a concise writing, it is useful to define 
some intermediate quantities; 

Ai(z) = A-z, (16) 
A2(z) = Az-K, (17) 



C/(z,7) = V(^'-l) (72-1) + ^7-1, (18) 
V{z, 7) = 2^(A2(z)-a2)(A2(^)_^2) + 2AA2(z)7 - 2(«:A2(z) + a^), (19) 

W±{z,j) = 2^{Al{z) - a2)(A2(7) - a^) ± 2XA,{z)j - 2{zA,{z) + a'), (20) 
e f \ fl27~\ ^1 (70 T 1)W^±(2,7) 

S±{z,-f,jo)^ JAi{z)-a^\n _^-,.„, . r, 21 

V (7T l)VF±(2:,7o) 

R(z..)- ^ In i^~"foP(^'"f) , VA2(z)-«2 (z^7o)T/(z,7) 

-K(^,7,7o) - ^===ln- — -H In- — -. (22) 

Vz^ - 1 (z - 7)[/(z,7o) z^-1 (z-7)F(z,7o) 



We can then write 

t 7;i?(7;,7,7o) + r/i?(-r/,7,7o) S'+(k,7,7o) S'_(-k,7,7o) 



(23) 



U A'(7,+rO 2A'(7,-l)(rz + l) 2A'(7j + i)(r, ~ 1) 

We can see on Fig. [21 how the speed of the ramjet tends toward the limit speed as time increases. 

It is interesting to look at the limit k — > 0. In this case, we have: 

11 — ' Y A^' ^ ^ 



p. ^ a/^, (25) 



A2(K)-a2^ Va'-A', (26) 



where a = e + A in the r.h.s. of these relations. Moreover, F/ and 7/ tend toward the same limit. In formula (|23p . 
one can see that, in the limit k — > 0, the coefficient of the functions R tends toward „ and the coefficient of the 

functions S± tends toward —7=—, where parameters 7/ and /3; are taken for k = 0. So, one can see that expressions 

1\ Oi Pi — 

for the function tjt^ obtained in this paper when the parameter k, vanishes and given by formula (19) in Ref. [7| are 
identical (the supplementary factor 1/2 can be absorbed in the In-functions). 
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FIG. 2: Reduced speed /3 of the ramjet as a function of the reduced time t/t, spent in the inertial frame of the interstellar 
medium, for the parameter u = 0.0071 [l[ and an initial reduced speed /3o ~ 0.1. The solid (dashed) curve corresponds to 
K = («: = 0.1). For each value of k, curves from bottom to top are respectively drawn for values 0.8, 0.5, 0.2 of the ratio X/uj. 



V. PROPER TIME 

Using the well known relation between the time and the proper time dr — dt/j [1, Q, Eq. (|15p simplifies and the 
proper time is given by the following integral 

1 r . r , ^. (28) 

t* Jo J JO (72 - 1) f ^72 - 1 + F{j) - y/j^~T] 

An analytical solution of this integral can be found with a procedure similar to the one used for the calculation of 
t/t^,. Using again the notations P^ - (P^ . a tedious calculation gives 

^ ^ R{li,l,lo) ~ ^(-ri,7,7o) _ '5'+(^:7:7o) 'S'-(-^:7:7o) ^20^ 

U X'hi+Ti) 2A'(7z-l)(r, + l) ^ 2A'(7, + l)(r,-l)- ^ ^ 

The link between the proper time on board the spaceship and the time spent in the inertial frame of the interstellar 
medium can then be computed (see Fig. [3]). 



Using relations (P^ - (|7f|) in formula (P^ . one can see that, in the limit k ^ 0, the coefficient of the functions R 
tends toward ^-^^ the coefficient of the functions S± tends toward ^^^fs ' ^^ere parameters 7; and Pi are taken 
for K — 0. Again, one can see that expressions for the function r/t^, obtained in this paper when the parameter k 
vanishes and given by formula (24) in Ref. f?] are identical (the supplementary factor 1/2 can be absorbed in the 
In- functions). 



VI. DISTANCE 

Using the relations 



Eq. (|15|) can be rewritten into the form 



^^%x=r — ^ ^. (31) 



* Jo J-ro 



t/t. 
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FIG. 3: Reduced time t/t* spent in the inertia! frame of the interstellar medium as a function of the reduced proper time 
r/t* on board the ramjet, for the parameter lo — 0.0071 1] and an initial reduced speed (3o ~ 0.1. The solid (dashed) curve 
corresponds to k — = 0.1). For each value of k, curves from bottom to top are respectively drawn for values 0.8, 0.5, 0.2 
of the ratio X/uj. 



Again, an analytical solution of this integral can be found with a procedure similar to the one used for the calculation 
of t/t^ and r/t*. New intermediate quantities must be defined: 



We can then write 



— =X(7)-X(7o) with 



X(7) = 



A'(7/+r/) 
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(32) 
(33) 

(34) 

(35) 

(36) 

(37) 

(38) 
(39) 



where 11 is the incomplete elliptic integral |10 |. The distance travelled by the ramjet in the interstellar medium as a 
function of the proper time indicated by the on board clocks can then be computed (see Fig.|4|). 



In the limit k ^ 0, the first term of formula p9p reduces to ^ argtanh where parameters 7; and Pi are taken 
for K = 0. Again, one can see that expressions for the function x/a;* obtained in this paper when the parameter k 
vanishes and given by formula (32) in Ref. are identical. 
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FIG. 4: Reduced distance a;/a;« travelled by the ramjet in the inertial frame of the interstellar medium as a function of the 
reduced proper time r/t, on board the ramjet, for the parameter lu — 0.0071 fl| and an initial reduced speed /3o ~ 0.1. The 
solid (dashed) curve corresponds to k = («: = 0.1). For each value of k, curves from bottom to top are respectively drawn for 
values 0.8, 0.5, 0.2 of the ratio A/a;. 



VII. NON-RELATIVISTIC LIMIT 



In the non-relativistic limit (/3 ^ 1), Eq. ^ reduces to 



where 

F(l) = e(2a + e). 

This number is clearly positive since e > and a < 1. Eq. ([40]) can be integrated to give 

X /3 - /3o 



with 



/3 = /3o exp ( v/e(2a + e) — 



(40) 

(41) 

(42) 
(43) 



These expressions, characteristic of an exponential motion 0,0], are valid provided t <C — ln/3o/ \/ e{2a + e). Let us 
remember that t does not differ significantly from t at non-relativistic speed. 



VIII. ASYMPTOTIC MOTION 



Near the limit speed, the term F{'j) is very small with respect to the term 7^ — 1, so we can write 



Equation ^ reduces then to 



- 1 + Fi^) - - 1 



if* 273 



(44) 



(45) 



It is worth noting that, in the case of an ideal ramjet {X = k = 0), ^(7) = e' and Eq. ((45|) is then the equation of 
motion of an uniformly accelerated spaceship with a constant proper acceleration ip^e' /2 |6|. 



If the ramjet is characterized by a mass M = 1000 t and an effective intake area A = 10^ km^, the characteristic 
acceleration is tp^, w 1500 m/s^ for an optimistic value of the particle density, let us say IQ-^/cm^ |l|. Within such 
conditions, the asymptotic proper acceleration of an ideal ramjet with e = 0.0071 is about 1 g. With an initial speed 
as low as 10 km/s, the spacecraft could then approach light velocity within a year. With a more realistic value, 1000 
times smaller, for the particle density, the performance of the ramjet would be reduced by the same factor. 

In the asymptotic regime, the reduced speed is very close the reduced limit speed. So we consider values of (3 such 
that < f3 < /3i, with (3i — (3s ^ Using Eq. ([8]), the factor 7^ = 7 (/3s) can then be approximately written 

7. = 7(A - [Pi - Ps)) ifM - A). (46) 

So, we have 

1 



Pi-Ps~—hi-ls). (47) 



The condition j3i — (3s (3i is then equivalent to 7; — 7s «C lf(3i < jf- This inequality is relevant since 7; has a finite 
value when A 7^ or k ^ 0. We assume that the ramjet is at position Xs, at time is, and at proper time Ts when it 
has a speed /3s- Using the approximation (HH), the integral for the position is given by 

~J.^ 4 A'(r,+7)(7i-7) A'(7,+ro 4 ^ ^ 

since 7^7/ under the integral. The corresponding solution is 

7 = 7/ - (7/ - 7s)exp I — (x-Xs)l. (49) 



One can treat the time exactly in the same way. Remarking that "f/\/j^ ~ 1 ~ 1/ (3i under the integral, the solution 
is given by 

7 = 7; - (7i -7s)exp I {t-ts)\. (50) 

This equation can be directly deduced from equation (j49p . since x — Xs ~ (Jic{t — ts) in this regime. Similar calculations 
for the proper time lead to 

/ X'(3ni{ji+Ti) ^ \ 
7 = 7i - (7i - 7s)exp I (r - Ts) 1 . (51) 

This equation can be directly deduced from equation (|50p . since t — ts ~ ^i(t — Ts) in this regime. 

Since (3s < (3 < (3i, Eq. I[i7|) holds also for (3s replaced by (3. Consequently, in Eqs. (|i5)) - ([5T|) . quantities 7, 7/, 7s 
outside the exponentials can be replaced respectively by /3, (3i, (3s- In the limit k — > 0, it can be checked that Eqs. (|49p . 
(fsn]) . (jSTj) tend respectively towards formulas (31), (17), (23) in Ref. 0] (note a misprint in Eq. (31): the factor (3i 
must be suppressed in the exponential). 

IX. FUNDAMENTAL INEQUALITIES 

Since 7 > 1, we have 

r=rdr=f'' " 
Jo Jo 

Moreover, since /? < 1, we have 



T = / dT= — < dt^t. (52) 
Jo Jo 7 Jo 



x= dx = (]cdt < cdt = ct. (53) 
Jo Jo Jo 

So, from these inequalities, we can conclude that functions ([23ll . p9|) and (|38p are characterized by 

r t 

-(7,70) < -(7,70), (54) 
X t 

— (7,70) < —(7,70)- (55) 



It is not evident to demonstrate these properties from the exphcit form of these functions, except in the case of a 
perfect antimatter ramjet (e = 1, X — k — 0) Nevertheless, we have checked that they are fulfilled numerically in 
any case. 



X. SUMMARY 



Formulas ([23|) . ([29|) and ((38|) form the complete set of parametric equations of motion for a Bussard ramjet with 
radiation and mass losses, as a function of its speed. It is then easy to compute, for instance, the distance traveled by 
the ramjet in the interstellar medium as a function of the proper time indicated by the on board clocks (see Fig.H]), 
or the link between this proper time and the time spent in the inertial frame of the interstellar medium (see Fig. [3]). 

With radiation and mass losses, the ramjet speed cannot be arbitrarily close to the speed of light. A limit speed, 
only reached asymptotically (see Figs. [T] and [21), lowers the performance of a Bussard ramjet as an interstellar spaceship 
or as a time machine for the exploration of the future (111] . 
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